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DUAL-RATE FINITE-SETTLING-TIME DISCRETE SYSTEMS 


by 

J. A. Gatlin 

Goddard Space Flight Center 


INTRODUCTION 

Finite-settling-time (FST) response for a single output/single input linear feedback control 
system can be obtained using the dual- rate, sampled-data algorithm 

n 

- u k+1 = b 0 (k) U k + ; 

i= 1 


u k = control level, kT^ t < (k + l) T, 
d. (k) = output error sampled at t = (k + i/n) t, 
b 0 (k), b. (k) = controller coefficients, i = 1 , 2, ■ • • , n; k = 0, 1 , 2, • • •. 

The performance thereby obtained is superior to that obtained using the FST control law for dis- 
crete systems discussed by Kalman and Bertram (Reference 1) and Lindorff (Reference 2). Their 
procedure is based upon the availability of all the state variables and a plant with a rational trans- 
fer function, and the design obtained suffers from high noise sensitivity and restricted dynamic 
range due to plant saturation. 

The control law serially accumulates n weighted, plant output error samples, 0. (k), 

i = 1, • * • , n, during each t second interval. Each sample is separated by T/n seconds, with the 

first sample taken at t = kT + (T/n) . At the end of the k th interval, that is, at t = (k + 1)T, the con- 
struction of a new control level is completed. Then, retaining only the value of the previous con- 
trol level, b 0 (k) u k , the control law construction cycle repeats. 

The synthesis method involves selecting the coefficients b. (k), i = 1, • • • , n, such that the 
system is brought to a stable null equilibrium in m + 1 control level updates, where m is the order 
of the plant. Also, the synthesis method minimizes 


f ■ !>■>' 

; = i 


(2) 


in order to minimize input noise transmission. 
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FORMULATION OF THE PROBLEM 


The differential equation describing the plant response for the system shown in Figure 1 can 
be written as 



x(t) - Gx(t) + gu( t) , 

u(t) = U k , kT < t < (k + 1) T , 


Figure 1— A dual-rate sampled data system. 


x(0) o x Q 


(3) 


where Gisamxm matrix, x( t ) is the plant m-state vector, x 0 defines the plant initial state, and g is 
the driving: vector which for this single input u(t) has all elements zero except the last which is the 
control gain constant K. Equation 3 has the solution 


x(t) - 0(t)x + <J>(t) 0)(“T) gu(T) dr , 


I' 


(4) 


where <I>(t) r e Gt is the plant transition matrix. Letting x k (-C) = x( t)| t=kT+ ^ T / n , Equation 3 yields 
the state vector at the error sample times during the k th interval as 


llr\ !ll\ & ,n 

(t) = ®\—) * k (°) + 4>(-T) gu(T) dr . 

^0 


(5) 


Selecting the first element of x( t ) as the plant output variable e yields the output samples 




( 6 ) 


where is the first row of the matrix a>, and 


-'0 


<J)(-r) g dr . 


(7) 


In what follows, <D(-tT/n) will be written as $('.). (fr/n) h(^T/n) , a scalar depending only on l, 
will be represented by &{Z) . For l = 1, 2, • • . , n, Equation 6 produces 


*k(l) 


(i) 


a(l) 


- 


x k <0) + 


K ( n )_ 


(n) _ 


_ a ( n )_ 


J k • 


(8) 
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Defining 


y k T - K (1) ’ * • ( n )] , n-vector , 

V T = ( 1) * • • </>* (n)] , n x m matrix , 

a T - [a(l) ••• a(n)] , n-vector , (9) 

reduces Equation 8 to 

y k = Vx k (0) + au k . (10) 

Likewise Equation 1 can be expressed in vector notation as 

" Vi = b o u k + bT yk ’ (11) 

where b T = [b, (k) • ■ • b n (k)] . Combining Equations 10 and 11 yields 

” u k+i = b T Vx k (0) + (b 0 +b T c) u k . (12) 

In Equation 12 the output samples, 0. (k), have been replaced by a function of the state of the 
plant at the beginning of the k th control interval, x k (0) , and control level, u k? applied during the 
interval. Let c, an m- vector, be defined as 


c T = b T V , 

c = V T b , (13) 


so that Equation 12 also has the form 


J k+i 


( 0 ) 


+ C 0 U k 


(14) 


where c 0 = b 0 + b T a. 

The control vector c T , c 0 can be synthesized from the weighting vector b Q , b T as follows. For 

n = m 


b = [VT]' 1 c . 

b o = c o ‘ bT ° • (15) 
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and for n > m + 1, b 0 can be left as a free parameter so that 



( 16 ) 


Thus for n = m + 1, 


b ( b o) = H-‘C m+1 , 

and for n > m + 1, the solution of Equation 16 that also minimizes 


(17) 


n 



i= 1 


b 0 PT ( b o) = H T [hH t ] _1 C m+1 • 
Also, if b Q is not left as a free parameter for n >m, 


(18) 


b oP t = V[V T V]-‘C , 

( b 0 ) OP T = C 0 " b OPT° ' ( 19 ) 

For either Equation 17 or 18, b 0 may be either arbitrarily selected, or selected so that F is also 
minimized with respect to b 0 . 

Minimization of the noise content factor, f, results in minimum noise transmission through the 
controller since for white noise input, 



(20) 


where o- 2 is the variance of the control level sequence and cr 2 is the variance of the output sample 
sequence. 

In this section the equivalence between Equations 1 and 14 has been developed. In the following 
section the unique values of c T , c 0 required for FST response will be determined. The values of 
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b 0 , b T used in Equation 1 are not unique but depend on the number and the format of the error 
samples. Equations 15, 17, 18, and 19 can be used to determine the values of b 0 , b T required to 
construct any given (m + 1) control vector c T , c 0 . 

DETERMINATION OF THE FST CONTROL VECTOR 

The change in the state of the plant during the k th interval is given by Equation 5 with l - n : 

X k (n) = * k+1 (°) = fl>(n)x k (0) + <D(n) h(n) u k . (21) 

Equations 14 and 21 can be combined to yield 


x k + l 

= A 

(0) 

" " 

A 


u k+l 


u k 


( 22 ) 


where 


<5(n) | $(n) h(n) 



[*u] 


(23) 


The matrix A is the controller plus plant, closed- loop, state transition matrix. To obtain FST re- 
sponse in (m + 1) steps for any set of initial conditions, the parameters c T , c 0 must be chosen so 
that 


A m+1 = 0 . (24) 

Direct solution of Equation 24 is quite tedious, but from the Caley-Hamilton theorem, Equa- 
tion 24 is satisfied if A, a (m + 1) x (m + 1) matrix, has all its eigenvalues equal to zero. The canon- 
ical form of such a matrix has zeros everywhere except the diagonal above the main diagonal which 
contains unity elements. Equation 24 is satisfied for 


IzI -A| 


z m + i 



0 , 


(25) 


if 


a. - 0, i “ 0, 1, • • ■ , m . 


(26) 
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There are (m + 1) equations involved in Equation 26 and (m + 1) controller parameters c T , c Q . in 
general, expansion of Equation 25 about the last row of i zi - Ai followed by the conditions of Equa- 
tion 26 yields 


Qc ~ q , 


c Q - Trace <I>(n) 


(27) 


The matrix Q and the vector q are systematically obtained from A as illustrated below for a 
third-order matrix: 


l z " a ll 

_ a i 2 

"“l j 



Z " a 22 

_a 2 l 


1 

C 2 

2 + c 0 ! 



Expanding about the last row yields 

0 = K a !2 -a i a 22 + “l Z ) c i + ( a i a 21 " a 2 a H +a 2 Z ) c 2 

+ ( a U a 22 “ a i 2 a 2 ! ‘ ( a il + a 2 2 ) Z +Z2 ) ( Z + C o) ' 

Setting a 2 = 0 produces c 0 = a u + a 22 , and setting a t = 0, a 0 = 0 yields 


(28) 


(29) 


0 


K 3 12 ' a i a 2 2 ) (‘ 


C 0 2 ' ( a il a 22 " a !2 3 2 1 ) 
; ( 3 11 a 22 _a !2 a 2l) C 0 _ 


(30) 


A SECOND-ORDER SYSTEM WITH REAL ROOTS 

Consider the system of Figure 1 with m = 2, a x =0, and 
a 7 = a. Given the parameters K, a, T,andn > 2, the controller 
parameters b 0 , b T are to be chosen to yield FST transient re- 
sponse and a minimum noise content factor, F. The fictitious 
controller parameters c T , c Q? illustrated in Figure 2, are 



Figure 2— A FST system with direct 
state variable feedback. 
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obtained first. Letting b = e aT results in 


4>(n) = 


>11 

^12 


1 

1 -b/a 

J21 

0 22 _ 


0 

b 


h(n) 


K 

a 


( 1 + aT - b" 1 )/a 

L ( b ' 1_1 ) - 


(3D 


so 


<D(n)h(n) - 


(aT + b - l)/a 
(1 - b) 


As defined in Equation 23 , the matrix A is 



1 

(i -b) 
a 

( aT + b - 1 ) , 

(aT) 2 KT 

A = 

0 

b 

( aT ) KT 


_c i 

“ C 2 

" C ° 


(32) 


using Equation 30 yields 


- O' 1 q 


aT 


(1 - b) KT 2 

1 - b - aTb 3 
(1 -b) 2 KT 


(33) 


with c 0 = 1 + b. It is sometimes useful to set c 1 = 1 ; that is, to use unity static position feedback. 
This allows the FST loop gain to be expressed as 


K 0 T 2 = (T^b) ' 


aT 


(34) 


and reduces the expression for c 2 to 


1 - b - aT b 3 
(1 -b)a 


( 35 ) 
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For n = 2, Equation 15 yields the weighting parameters b QJ b 1} b 2 from the calculated ficti- 
tious controller parameters c 0 , c 1 , c 2 . Using 


V T 


1 1 

1 - b 0, 5 1-b 

a a 


(36) 


and 


KT 2 aT/2 + b 0, 5 - 1 
( aT ) 2 aT + b - 1 


(37) 


consider the case for aT^ 0. Since e * = l - x + x 2 /2! - x 3 /3! + ■ • • , the general expressions reduce 
to 


< 0 T 2 " at/faT- \ (aT) 2 + 


'aT- j ( aT ) 2 + **’ " a T + 3(aT) 2 


aT - \ (aT) : 


5T 

2 


with c 0 = 2. Equation 15 now is 


- 1 

1' 

- 1 

- i - 


--3" 

_T/2 

T- 


- 5T/2 _ 


_ 4_ 


b 0 = 2 - [-3 4] 


1/8 
.1/2 J 


= 3/8 


(38) 


(39) 


For n = 2 the noise content factor is 

F = (4) 2 + ( 3) 2 = 25 . 

Now consider the use of n = 3. The fictitious controller parameters, for KT 2 = 1, are 


(40) 



" *i 1 


1 

^m+l 

C 2 

= 

5T/2 


o 

I 

cr 

o 

1 


_2- b o 
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and as defined in Equation 16, 


Equation 17 yields 


H = 


1 

T/3 

1/18 


1 

2T/3 

4/18 


1 

3T/3 

9/18 


9 

b i 4 9 b o ■ 

b 2 = - 9 + 18 b 0 

b 3 = T ' 9b 0 • 


Calculating the noise content factor and dF(b 0 )/clb 0 


= 0 produces 



resulting in (b 0 ) OPT = 0.5185, F MIN = 15.5, and 


b 


OPT 


-2.417 
0. 333 
3.083 


(42) 


(43) 


(44) 


(45) 


A unique feature of the general control law expressed by Equation 1 is the incorporation of 
memory of the previous control level. If b 0 - 0 is used, FST response is obtained at the expense 
of F = 146. The low noise transmission feature of the disclosed FST control law is one of its most 
attractive features. 

For n > m + 1 the matrix H is (m + 1) x n. Equation 18 yields, for n = 4, 


“ll - 32 b~ 


”2 

-21 + 32 b Q 


-0. 5 

-15 +32b 0 


1 

29 - 32 b 0 i 


2. 5 


(46) 


with ( b o ) opt = 0.5938 andF MIN = 11.5. 
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For n > m + 1 it is sometimes impractical to use n distinct coefficients. Using only four, Equa- 
tion 1 takes the form 


J k+i 


- b n u,. 


n 7 4 n/2 3n/4 n 

+ b, 22 (k) + b 2 0i (k) + b 3 9 1 (k) + b4 H &i (k) 


i=n/4+l 


i=n/2 + l 


i=3n/4+l 


(47) 


With the error samples thus grouped in 4 sets of n/4 samples, the matrix H is 


H 



_n 

4 


n 

T 


n 


n + 4 3n + 4 5n + 4 

32 T 32 T 32 T 


7n + 4 
32 T 


n 2 + 6n + 8 
384n 


7n 2 + 18n +8 
38 4n 


19n 2 + 30n +8 
38 4n 


37n 2 + 42n + 8 
384n “ 


and for large n, H becomes 


H 



1 

T/8 

1/96 


1 1 

3T/8 5T/8 

7/96 19/96 


1 

7 T/8 
37/96 


resulting in 



-2.150 
-0. 550 
1.050 
2.650 


with (b 0 ) opT = 0.8333 and f min = (13.05) (4/n). 
As n - ^Equations 47 and 49 become 


(48) 


(49) 


(50) 


b o u k + 


/•T/4 

h 

•' n 


6 { k. t) dt + b. 


r 

J T /4 


0 ( k, t) dt + b. 


/• 3 T • ' 4 
^T/2 


#(k , t ) dt + b 4 


r 

■' 3 T/4 


*(k. t)dt ; (51) 


H 



1 

3T/8 


1 

5T/8 


1 

7 T/8 


1/96 7/96 19/96 37/96 


(52) 
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with (t/ 4) b opx = (n/4)b OPT . The calculation of f min , however, must now include the description of 
the noise source power spectrum. As an example, consider a band-limited noise generator with 
uniform power density, 


p d = 


2 /f - f < f < f 

o/ I o ’ r o 1 I 0 ’ 


(53) 


passed through a linear, low-pass filter with the Laplace transfer function 


jfK 

G(s) = \ « 2 tt f 0 . 


(54) 


The resulting noise source is described by the variance 


1 r p <* i , ./ , 

= ~xn I T G(s)G(-s)ds = "4 (^o/ f 0 ) 

J J- \ 00 


(55) 


and the auto -correlation function, 


R Nn( T ) = (56) 

For samples taken every T/n < 1 A seconds, it is no longer appropriate to assume statistical inde- 
pendence of the samples; however, if the noise source is integrated for T/4 > 1 A seconds, samples 
taken every t/ 4 seconds can be assumed independent with 


/•T/4 

■ J (5-0 

n 


R nn ( T ) dT ’ 


(57) 


which yields 


T 

4x: 


(58) 


CT 2 CT 2 


[ b 4 2 +b 3 2 + b 2 2 +b l 2 ] ’ 

e I e 

Therefore, as for Equation 50, f min = (13.05) (4/\t). 


( 59 ) 
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The preceding cases show how F is reduced for n > m; Table 1* is given in summary. 


Table 1 

Effect of Sample Averaging. 


n 

(ko)oPT 


AT/n 

2 

0.3750 

5.000 

1024 

3 

0.5185 

3.937 

683 

4 

0.5938 

3.391 

512 

32 

0.8022 

1.267 

64 

128 

0.8255 

0.637 

16 

512 

0.8314 

0.319 

4 

CO 

0.8333 

0.160 

0 


DESIGN EXAMPLES 

Consider the system of Figure 1 with m = 2 , a 1 = 0, a 2 = 0. The resulting K/s 2 plant model 
is used as the basic model for many spacecraft single-axis attitude control designs. Using the pro- 
cedure given in the preceding sections, the matrix A is 



1 

1 

0.5KT 2 

A = 

0 

1 

KT 2 


-1 

-2.5 

( 

cr 

o 

1 

to 

I 

cr 

o 

¥ 

to 



(60) 


with the c T , c 0 parameters chosen to yield FST response for KT 2 = 1. A root locus plot of a(kt 2 ), 
given in Figure 3, illustrates the relative stability of this FST design. Note that for KT 2 = 0 there 
is a z-domain pole at z = - b 0 ; therefore if b 0 < 1, system stability for K max < K < 0 is obtained. 


For n = 2 the control algorithm is 


- u k+1 = 40 2 (k) - 3B X (k) + (3/8) U k 


(61) 


*For noise with R NN 


e Kr , AT = 2048. 
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with F = 25; for n = 4, F min = 11.5 for (b 0 ) OPT = 0.5938 
and 

b 0 T PT = [~2 -0.5 1 2.5] . (62) 

Figure 4 shows how the elements of b 0PT vary for n = 10 
as a function of b 05 and Figure 5 shows the variation of F 
with respect to b 0 . 

If only 3 distinct coefficients, b 2 , b 1? b 0 , are used 
with n = 10, the control algorithm is 


-"k+i 


= (4.4/5)^V <k) -(3.4/5) &. (k) +0 . 735u k ■ (63) 



Figure 3— FST roof locus for b 0 = 0.8333. 




Flgure4— The elemenfs of bop T for n = 10. Figure 5— F versus b 0 for n - 10. 

with F = [(4.4) 2 + (3.4) 2 ]/5 = 6.18. Thus, for n = 10, use of sample averaging to reduce the num- 
ber of required distinct controller parameters increases the noise factor F by 6.18/4.71 or 31%. 
For large n, nb 2 ^2 = 4.5, nb 1 ^2 = - 3.5, and F = 85/n. 


Consider now the use of pulse- width- modulation during the T second interval so that re- 
action jet attitude control torquers can be used. The jet on-time AT k is scaled to apply the 
same control impulse as the preceding amplitude -modulation design; that is, 


AT u 

k max 


u kT 


AT k < T . 


(64) 
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The matrix A, now non-linear, becomes 


1 

1 i ( 

i 

1 - 0.5y k ) KT 2 



0 

■ 

i : 

i 

• 

KT 2 

0 < y k < 1 . 

(65) 

3 

-c 2 : 

_c o _ 




y k = AT k/ T = K|/ U m a* ' (66) 

Applying the linear FST design methods, assuming y k is constant, yields finite settling time 
response as y k -constant with c 0 = 2, = 1, and c 2 = 2 + 0.5 y k . Using Equation 19 and 3 dis- 

tinct controller coefficients yields 


/ n \ 

~-( 2 - 5+ yJ 


_ 3.5 + y k _ 


b 0 = 2 - b T ec(r k ) . 


The elements of a are piecewise continuous functions of y k ; for example, 


a 


12 


3 - 6y k + 4y k 2 
9-6r k 


0 < y k < 0.5 . 


(67) 


( 68 ) 


Also for y k = 0, b 0 = 0; and for y k = 1, b and b 0 are the same as used with amplitude modulation. 

The matrix A can be used to examine the non-linear response of this FST design. Defining u k 
as the calculated control level and retaining u k as the control level applied and remembered allows 
minimum jet on-time and jet full-on conditions to be handled by the logic 


U k = 0 , 


0 < u. < U . 

k min 


U . < u. I < U 

mm Ik max 


u. - U 

k max 


U < U. 

max k | 


(69) 


For large initial rates the response for KT 2 = 1, u max = 1, U min = 0.05 is given in Tables 2 and 3. 
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Table 2 


Table 3 


Saturated Transient Response. 


k 

*k 

T0 k 

u k 

I 

u k 

0 

0 

1 

0 

0 

1 

1 

1 

-1 

-2 

2 

1.5 

0 

-1 

-1.5 

3 

1 

-1 

0.5 

0.5 

4 

0.375 

-0.5 

0.25 

0.25 

5 

0.093 

-0.25 

0.187 

0.187 

6 

0.013 

-0.063 

0.054 

0.054 

7 

0.004 

-0.008 

0 

0.003 

8 

-0.004 

-0.008 

o 

0.011 

9 

-0.012 

-0.008 

0 

0.019 

10 

-0.02 

-0.008 

0 

0.026 


Unsaturated Transient Response. 


k 

*k 

T0 k 

u k 

u k 

0 

0 

0.25 

0 

0 

1 

0.25 

0.25 

-0.5 

-0.5 

2 

0.125 

-0.25 

0.187 

0.187 

3 

0.044 

-0.063 

0 

0.023 

4 

-0.018 

-0.063 

0.08 

0.08 

5 

-0.004 

0.017 

0 

-0.016 

6 

0.014 

0.017 

0 

-0.032 

7 

0.031 

0.017 

0 

-0.05 

8 

0.049 

0.017 

-0.068’ 

-0.068 

9 

0.002 

-0.05 

0 

0.048 

10 

i -0.048 

! -0.05 

0.096 

0.096 


SIMULATION RESULTS 


Construction of a FST system to control a second-order plant was accomplished using hybrid 
simulation. The continuous plant, G(s) - K/s 2 , was modeled on the analog computer, and digital 
computation generated the control law of Equation 1 and the output error 6 = 6 0 - where 0 Q is 
the plant output and 0 R is a staircase reference function that changes every T seconds. 

The control law update interval, T, is 11 seconds, of which 10 seconds is used for output sam- 
pling and 1 second for the computational and conversion delays; that is, Equation 1 has the form 


U k + 1 = b o U k + b 


n/4 n/2 3n/4 

! + b 2 (1<) + b3 2l (k> 


i=n/4 + l 


+ b. 


n 1 1 /At 

22 e i (k> + b 5 21 (k) 


i=3n/4+l 


(VO) 


with b 5 = 0 and n = (T-TJ/&T = 10/AT . Output error sampling periods, AT, of 10/512, 10/128, and 
10/32 seconds are provided; an option allowing the plant output to be integrated for 10/4 seconds 
on the analog computer and sampled every 10/4 seconds is included. 


15 





Illlllllllllllll 


Since the data sampling interval is now only (t -T d ) , Equation 49 becomes 


_ c o b o_J 


■ (t) 


(T-TJ 

8 

(t-T,) 2 /t 2 

96 


3(T-T d ) 

8 

7(T -T d )^T 2 
96 


S(T-T d ) 


96 


7(T-T d ) 


19(t - T d ) 2 /t 2 37 (t - T d ) /t 2 


96 


( 71 ) 


and letting d = T d/ /(T - T d ) , Equation 71 reduces to 


C 1 


1111 







b. 

(1 + d)c 2 T 

■ (i) 

1/8 3/8 5/8 7/8 


2 

b 3 

_0 +d ) S ( c o ' b o) 


! 1/96 7/69 19/96 37/96 


_ b 4_ 


(72) 


As shown in Reference 3, the relative stability as described by the root locus plot (Figure 3) 
depends only on the plant gain, K, and the value of b 0 . If b 0 is kept at the value of (b 0 ) OPT for d = 0, 
the addition of the delay affects only the value of F. The inverse of Equation 72 yields 
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for b 0 = 0.8333. For d = 1/10, 
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0. 190 
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(73) 


(74) 


with F = (20.13) (4/n) which corresponds to an increase in F of 54%. 

As developed previously (Equations 34 and 38) the FST loop gain requires a plant constant, K 0 , 
such that 
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If K / K 0 , a non-finite but rapidly convergent sequence will characterize the transient response. 
Figure 6 presents this effect for an initial condition response, d Q £ 0, 0 o = 0; Figures 7, 8, 9 show 
this effect when the controller is tracking an input "staircase" sinusoid. 


The coefficient set of Equation 74 is exact 
for n - oo and i S a good approximate set for large 
n. Figure 10 shows the effect of the approxima- 
tion for n = 32, n = 128, and n = 512; Fig- 
ure 11 compares the sampling format for finite 
and infinite n. 

The steady-state performance when output 
noise is present shows, in Figure 12, the re- 
duction of the controller noise sensitivity as the 
error sampling rate is increased. The runs in 
Figure 10 correspond to the values of n = 32, 
128, 512, and 5 given in Table 1 of the preceding 
section. Figure 13 presents a transient re- 
sponse with noise present for n = 512. 


TIME 



Figure 7— Tracking response for K = 0.8 K 0 . 
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(a) K=0.8 K q (b) K=K q (c) K = 1 . 1 K q 

Figure 6 — Initial rate transient- response. 




Figure 8— Tracking response for K = K Q . 
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min i 


TIME - 



* RECORDER SATURATION 


Figure 9— Tracking response for K - 1.1 K 0 . 


CONCLUSIONS 

This paper has extended the work done in 
Reference 3 so that the dual-rate FST control 
law can be applied to n th order single-input, 
single-output systems. A means of handling 
computation and conversion delays, additional 
noise analysis, and additional hybrid simulation 
work is also documented. 

The general control law formulation given 
in Equation 1 allows for the controller coeffi- 
cients b 0 (k), b. (k) to be a function of the con- 
trol interval index, k. In one of the design ex- 
amples presented in this report, the coefficients 
required are constant. The need for varying 
coefficients arises when the matrix A (Equa- 
tion 23) has elements that change from control 
interval to control interval. This requires that 
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( a ) n = 32 


(b) n = 128 


( c ) n = 512 


Figure 10 — The effect of the large n assumption. 
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Figure 11— Error sampling format. 
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Figure 13 — Transienf performance with noise. 


*** N + q -<no!*-*) 


the control vector c T , c 0 also be a function of k 
Figure 12— Steady state performance with noise. if A is to approach an FST matrix as its ele- 

ments approach a constant value. This approach 
has been used to apply the FST control law to pulse-width-modulation control of a second-order 
plant. The elements of A can also change due to changes in the equilibrium position of a non-linear 
plant. Thus, using the non-linear model, the coefficients b 0 (k) , b i (k) can be adjusted to Tr track 
the operating point' T and maintain FST response about the operating point. 


In summary the key features of the new FST control law are: 

(1) The dual error sampling/control level update rates allow the error sampling rate to be 
chosen to meet noise sensitivity requirements and the control update rate chosen to avoid 
plant saturation. 

(2) The memory of the previous control level greatly reduces the need for the controller coef- 
ficients to produce derived error rate information. 


(3) Only the output state need be directly sensed; all the state variables appear combined in 
the generated control level, but their isolation is not required. 
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(4) The control technique possesses null stability and reasonable tolerance to errors in the 
plant model. 

(5) The finite memory of the control law allows the controller weighting coefficients to be 
readily changed during each control interval. 

(6) A multiplicity of error samples taken while the control level remains constant provides a 
source of data for plant and disturbance on-line modeling. 
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